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ABSTRACT 

We use a suite of cosmological N-body simulations to study the properties of substructure halos 
(subhalos) in galaxy-sized cold dark matter halos. We extend prior work on the subject by considering 
the whole population of subhalos physically associated with the main system. These are defined as 
subhalos that have at some time in the past been within the virial radius of the halo's main progenitor 
and that have survived as self-bound entities to z = 0. We find that this population extends beyond 
three times the virial radius, and contains objects on extreme orbits, including a few with velocities 
approaching the nominal escape speed from the system. We trace the origin of these unorthodox orbits 
to the tidal dissociation of bound groups of subhalos, which results in the ejection of some subhalos 
along tidal streams. Ejected subhalos are primarily low-mass systems, leading to mass-dependent 
biases in their spatial distribution and kinematics: the lower the subhalo mass at accretion time, 
the less centrally concentrated and kinematically hotter their descendant population. The bias is 
strongest amongst the most massive subhalos, but disappears at the low-mass end: below a certain 
mass, subhalos behave just like test particles in the potential of the main halo. Overall, our findings 
imply that subhalos identified within the virial radius represent a rather incomplete census of the 
substructure physically related to a halo: only about one half of all associated subhalos are found 
today within the virial radius of a halo, and many relatively isolated halos may have actually been 
ejected in the past from more massive systems. These results may explain the age dependence of 
the clustering of low-mass halos reported recently by Gao et al, and has further implications for (i) 
the interpretation of the structural parameters and assembly histories of halos neighboring massive 
systems; (ii) the existence of low-mass dynamical outliers, such as Leo I and And XII in the Local 
Group; and (iii) the presence of evidence for evolutionary effects, such as tidal truncation or ram- 
pressure stripping, well outside the traditional virial boundary of a galaxy system. 
Subject headings: cosmology: dark matter - methods: N-body simulations - galaxies: kinematics and 
dynamics - galaxies: halos 



1. INTRODUCTION 

In the current paradigm of structure formation, the 
concordance ACDM scenario, the dark matter halos 
that host galaxy systems are assembled hierarchically, 
through the merger and accretion of smaller subunits. 
One relic of this process is the presence of substructure, 
which consists of the self-bound cores of accreted subsys- 
tems that have so far escaped full disruption in the tidal 
field of the main halo (Klypin et al 1999; Moore et al 
1999). 

Although substructure halos (referred to hereafter as 
"subhalos" , for short) typically make up only a small 
fraction (5 to 10%) of the total mass of the system, they 
chart the innermost regions of accreted subsystems, and 
are thus appealing tracers of the location and kinematics 
of the galaxies that subhalos may have hosted. Substruc- 
ture is thus a valuable tool for studying galaxies embed- 
ded in the potential of a much larger system, such as 
satellite galaxies orbiting a primary, or individual galax- 
ies orbiting within a group or cluster of galaxies. 

This realization has prompted a number of studies over 
the past few years, both analytical and numerical, aimed 
at characterizing the main properties of subhalos, such as 



their mass function, spatial distribution, and kinematics 
(e.g. Ghigna et al 1998, 1999; Moore et al 1999; Taylor 
& Babul 2005a,b; Benson 2005; Gao et al 2004; Diemand 
et al 2007a,b). 

Consensus has been slowly but steadily emerging on 
these issues. For example, the moss function of subhalos 
has been found to be rather steep, dN/dM oc M^]^^ or 
steeper, implying that the subhalo population is domi- 
nated in number by low-mass systems but that most of 
the substructure mass resides with the few most massive 
subhalos (Springel et al 2001, Helmi, White & Springel 
2002, Gao et al 2004). Confirmation of this comes from 
the fact that the total fraction of mass in subhalos is 
rather low (typically below 10%) even in the highest res- 
olution simulations published so far (although see Die- 
mand et al 2007a,b for a differing view). 

Subhalos have also been found to be spatially biased 
relative to the smooth dark matter halo where they are 
embedded, avoiding in general the innermost regions. 
Furthermore, the number density profile of the subhalo 
population also differs markedly from that of galaxies 
in clusters, and possibly from the radial distribution of 
luminous satellites around the Milky Way (Kravstov et 
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al 2004; Willman et al 2004; Madau et al 2008). This 
precludes identifying directly the population of "surviv- 
ing" subhalos with galaxies in clusters, and highlights 
the need for either more sophisticated numerical model- 
ing techniques, or for pairing up the N-body results with 
semi-analytic modeling in order to trace more faithfully 
the galaxy population (Springel et al 2001, De Lucia et 
al 2004, Gao et al 2004, Croton et al 2006). 

One intriguing result of all these studies has been the 
remarkably weak dependence of the properties of sub- 
structure on subhalo mass. Gao et al (2004) and Die- 
mand, Moore & Stadel (2004), for example, find that 
the radial distribution of subhalos is largely independent 
of their self-bound mass. This is surprising given the 
strong mass dependence expected for the processes that 
dictate the evolution of subhalos within the main halo, 
such as dynamical friction and tidal stripping. Although 
efficient mixing within the potential of the main halo is a 
possibility, an alternative explanation has been advanced 
by Kravtsov, Gnedin & Klypin (2004). 

These authors argue that the present-day mass of a 
subhalo may be a poor indicator of the original mass of 
the system, which may have been substantially larger at 
the time of accretion, and used this idea to motivate how 
the faintest dwarf companions of the Milky Way were 
able to build up a sizable stellar mass through several 
episodes of star formation despite their shallow present- 
day potential wells. The same idea was also adopted 
by Libeskind et al (2007) as a possible reason for the 
peculiar spatial alignment of satellites around the Milky 
Way (Lynden-BeU 1976, 1982; Kunkel & Demers 1976; 
Kroupa, Thies & Boily 2005). 

We revisit here these issues with the aid of a suite of 
high-resolution N-body simulations of galaxy-sized halos. 
We extend prior work by carefully tracking the orbits of 
surviving subhalos back in time. This allows us to se- 
lect a complete set of subhalos physically associated with 
the main halo, rather than only the ones that happen to 
be within the virial radius at a particular time. As we 
discuss below, a large fraction of the associated subhalo 
population are on unorthodox orbits that take them well 
beyond the virial radius, a result with important implica- 
tions for studies of satellite galaxies and of halos clustered 
around massive systems. 

The plan of this paper is as follows. We introduce 
briefly the numerical simulations in § [2] and describe our 
subhalo detection algorithm and tracking method in § [3l 
Our main results are presented in § |4l we begin by ex- 
ploring the subhalo spatial distribution and kinematics, 
as well as their dependence on mass, and discuss the con- 
sequences of our findings for the subhalo mass function. 
We end with a brief summary and discussion of possible 
implications and future work in § [5j 

2. THE NUMERICAL SIMULATIONS 

2.1. The Cosmological Model 

All simulations reported here adopt the concordance 
ACDM model, with parameters chosen to match the 
combined analysis of the first-year WMAP data release 
(Spergel et al 2003) and the 2dF Galaxy Redshift Sur- 
vey (Colless et al 2001). The chosen cosmological pa- 
rameters are Vt^ = r^dm -I- f^b = 0.25, fib — 0.045, 
h = 0.73, r^A = 0.75, n = 1, and CTg = 0.9. Here Q. de- 
notes the present-day contribution of each component to 



the matter-energy density of the Universe, expressed in 
units of the critical density for closure, p„it = 3H^/8nG; 
n is the spectral index of the primordial density fluc- 
tuations, and (Ts is the rms linear mass fluctuation in 
spheres of radius 8 h~^ Mpc at z = 0. Bubble's "con- 
stant" is given by H{z) and parameterized at 2; = by 
H{z ^0) = Ho = 100 /i km s^^ Mpc^i. 

2.2. The Runs 

Our analysis is based on a suite of 5 high-resolution 
simulations of the formation of galaxy-sized ACDM ha- 
los. The simulations target halos of virial mass"'^, M200 ~ 
10^^ Mq, and have at z = between 3 and 5 mil- 
lion particles within the virial radius, 7'2oo- Each halo 
was selected at random from a list of candidates com- 
piled from a cosmological N-body simulation of a large 
(100 ft.^^ Mpc) periodic box and resimulated individu- 
ally at higher resolution using the technique described 
in detail by Power et al (2003). We imposed a mild iso- 
lation criterion (that no neighbors with mass exceeding 
5 X 10"/i-1Mq be found within Ih'^ Mpc at z = 0) 
in order to exclude systems formed in the periphery of 
much larger groups or clusters. 

The simulations were run with Gadget 2, a massively- 
parallel cosmological N-body code (Springel 2005). Par- 
ticle pairwise interactions were softened using the "op- 
timal" gravitational softening length scale suggested by 
Power et al (2003); i.e., a spline lengthscale = lAec ~ 
4 r2oo/\/ A'200 , kept fixed in comoving coordinates. Nu- 
merical details of each run are listed in Table [TJ 

3. THE ANALYSIS 

3.1. Substructure Finding 

We use SUBFIND (Springel et al 2001) in order 
to identify self-bound structures in N-body simula- 
tions. SUBFIND finds substructure within friends-of- 
friends (FOP; Davis et al 1985) associations by locating 
overdense regions within each FOF halo and identifying 
the bound subset of particles associated with each over- 
density. SUBFIND also works recursively and its output 
readily identifies "subhalos within subhalos" , thus char- 
acterizing fully the various levels of the hierarchy of sub- 
structure present within a given FOF halo. We retain 
for our catalogue all SUBFIND subhalos with more than 
20 particles. 

The main output of SUBFIND is a list of subhalos within 
each FOF halo, together with their structural properties. 
For the purposes of this paper, we shall focus on: (i) the 
subhalo self-bound mass, Mgnb', (ii) the peak of its circu- 
lar velocity profile (characterized by rmax and Knax); and 
(iii) the position of the subhalo center, which we iden- 
tify with the particle with minimum gravitational poten- 
tial energy. We have run SUBFIND on all 100 snapshots 
(equally spaced in scale factor, a) of each of our runs, 

^ We define the virial mass of a halo, M200, as that con- 
tained within a sphere of mean density 200 X Pcrit- The virial 
mass defines implicitly the virial radius, r"2oo, and virial velocity, 
V200 = y/GM2oo/r200t of 3, halo, respectively. We note that other 
definitions of "virial radius" have been used in the literature; the 
most popular of the alternatives adopts a density contrast (relative 
to critical) of A fs; 178 f^JJ;''^ ~ 100 (for our adopted cosmological 
parameters, see Eke et al 1996). We shall refer to these alternative 
choices, where appropriate, with a subscript indicating the value 
of A; i.e., rioo is the virial radius obtained assuming A = 100. 
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and are therefore able to track in detail the evolution of 
individual subsystems and their particle members. 

3.2. Substructure Tracking 

Our analysis focuses on all surviving subhalos at z = 
and relies heavily on tracking accurately their accretion 
history. To this aim, we trace each subhalo backwards 
in time by identifying the central particle at z = and 
searching for the group it belongs to in the immediately 
preceding snapshot. A new central particle is then se- 
lected and the procedure is iterated backwards in time 
until z = 9, the earliest time we consider in the analysis. 

This procedure leads in general to a well-defined evolu- 
tionary track for each subhalo identified at z = 0. When 
no subhalo is found to contain a subhalo's central particle 
in the immediately preceding snapshot, the search is con- 
tinued at earlier times until either a progenitor subhalo 
is found or z = 9 is reached. This is necessary because a 
subhalo may temporarily disappear from the catalogue, 
typically at times when it falls below the minimum parti- 
cle number or else when it is passing close to the center of 
a more massive system and its density contrast is too low 
to be recognized by SUBFIND. Our procedure overcomes 
this difficulty and in most cases recovers the subhalo at 
an earlier time. We note that these complications are a 
fairly common occurrence in the analysis procedure, and 
we have gone to great lengths to make sure that these 
instances are properly identified and dealt with when 
constructing our subhalo catalogue and their accretion 
histories. 

The tracking procedure described above defines a 
unique trajectory for each subhalo identified at z = 0. 
This trajectory may be used to verify whether a subhalo 
has, at any time in the past, been accreted within the 
(evolving) virial radius of the main halo. If this is the 
case, we record the time it first crosses r2oo{z) as the 
"accretion redshift", Zacc, and label the subhalo as asso- 
ciated with the main system. Analogously, we identify a 
set of associated dark matter particles by compiling a list 
of all particles that were at some time within the virial 
radius of the main halo but are not attached to any sub- 
structure at z = 0. On the other hand, halos that have 
never been inside the virial radius of the main halo will 
be referred to as "field" or "infalling" halos. 

Using the subhalo trajectories, we compute and record 
a few further quantities of interest for each subhalo; 
namely, 

• its "turnaround" distance, rta, defined as the max- 
imum separation between a subhalo and the center 
of the main progenitor before z = Zacc (for associ- 
ated subhalos) or before z ^ (for field subhalos); 

• the structural properties of associated subhalos at 
z = and at accretion time, such as mass and peak 
circular velocity; 

• an apocentric distance, Tapo, defined as the apoc- 
enter of its orbit computed at z = using the sub- 
halo's instantaneous kinetic energy and orbital an- 
gular momentum, together with the potential of 
the main halo^ 

^ We have checked that our results are insensitive to the triaxial 



Subhalo quantities measured at accretion time will be 
referred to by using the sub/superscript "acc" ; for exam- 
ple, V^ax refers to the peak circular velocity of a subhalo 
at z = Zacc- Quantities quoted without superscript are 
assumed to be measured at z = unless otherwise spec- 
ified; e.g., Knax = Knax(2; = 0). 

4. RESULTS 

The basic properties of our simulated halos at z = are 
presented in Table [H Here, eG(= ^s/1-4) is the Gadget 
gravitational softening input parameter, and M200: ^200 j 
and A^200 are, respectively, the halo virial mass, radius, 
and number of particles within r2oo- Table [T] also lists 
the peak of the circular velocity of the main halo, Knax, 
and its location, r„iax; the total number of "associated" 
subhalos; as well as the number of those found within 
various characteristic radii at z = 0. 

4.1. Subhalos beyond the virial radius 

One surprise in Table [T] is that the number of "asso- 
ciated" subhalos exceeds by about a factor of ^ 2 the 
total number of subhalos identified within r2oo ■ This re- 
sult is also illustrated in Fig.[Tl where we show, at z = 0, 
the distance from the main halo center vs radial velocity 
for all subhalos identified in our simulations. Distances 
and velocities have been scaled to the virial quantities of 
each primary halo. Colored dots are used to denote "as- 
sociated" subhalos, black symbols for "field" halos. Dif- 
ferent colors correspond to different subhalo masses, as 
measured by the peak circular velocity at accretion time 
(in units of the present-day primary halo virial velocity, 
V200): red is used for subhalos with V^ax — 0.72 V200, 
blue for those with F^^^ax — 0.038 V200, green for the rest. 

Note that the distribution of associated subhalos ex- 
tends well past ^ 3r2oo; indeed, a few associated subsys- 
tems are found at r ~ 4 r2oo moving outwards with radial 
velocity of order Vr ~ V2oo- A careful search shows that 
there are actually several associated subhalos presently 
at distances larger than ~ 5r2oo- 

This result is unexpected in simple formation scenar- 
ios, such as the spherical secondary infall model (SSIM, 
for short). SSIM identifies at any time three distinct 
regions around a halo: (i) an inner "virialized" region 
where accreted mass shells have had time to cross their 
orbital paths; (ii) a surrounding "infall" region, where 
shells are still on first approach and have not yet crossed; 
and (iii) a still expanding outer envelope beyond the 
current turnaround radius (Fillmore & Goldreich 1984; 
Bertschinger 1985, White et al 1993, Navarro & White 
1993). 

One of the premises of the secondary infall model is 
that the energy of a mass shell accreted into the halo is 
gradually reduced after its first pericentric passage until 
it reaches equilibrium. During this process, the apocen- 
tric distance of the shell is constantly reduced; for ex- 
ample, taking as a guide the SSIM self-similar solutions 
of Bertschinger (1985), the second apocenter of an ac- 
creted shell (the first would be its "turnaround" radius) 
is roughly 90% of its turnaround distance, rta, and the 

nature of the halo by recomputing rapo using the potential along 
each of the principal axes of the halo's mass distribution rather 
than the spherical average. This leads to typical variations of loss 
then ~ 20% in rapo- 
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Fig. 2. — Turnaround radius versus apocentric distance at 2: = 
for all associated subhalos in our simulations. The turnaround dis- 
tance is the maximum distance from the main halo before accretion. 
Subhalos on "traditional" orbits are expected to have rapo < fta 
and, thus, to be to the left of the 1:1 curve in this plot. Subha- 
los near the 1:1 line have rapo ~ rta and are therefore on orbits 
which have not been decelerated substantially since turnaround. 
Subhalos with rapo > rta are on unorthodox orbits and they have 
gained orbital energy during or after accretion. The blue symbols 
in the panel highlight subhalos on extreme orbits, that will take 
them more than ~ 2.5 times farther than their turnaround radius. 
The fraction of associated subhalos and associated dark matter 
particles in each region of the plot is given in the legends. 



riodic orbit with rapo ^ 0.8 rta. Thus, according to the 
SSIM, few, if any, associated subhalos are expected to 
populate the region outside 0.8 rta ~ 1.6r2oo- This 
is clearly at odds with the results shown in Fig. [T] and 
Table [TJ more than a quarter of all associated subhalos 
are found beyond 1.6 r2oo at z = 0! 



Fig. 1. — Upper panel: Radial velocity versus distance to the 
main halo center for all subhalos within 5 X r2oo in our simulations. 
Velocities and distances are normalized to the virial velocity, V200 1 
and virial radius, r200i of each host. All "associated" halos are 
shown in color, subhalos on first infall are shown in black. Differ- 
ent colors are used according to the peak circular velocity of the 
subhalo at the time of accretion. Blue denotes the quartile with 
smallest V^'^, red those with largest V^'^. Green denotes the 
rest of the associated subhalo population. Upward vertical arrows 
of matching color indicate the half-number radius for the various 
subhalo populations. A shorter black arrow marks the half-number 
radius for "associated" dark matter particles. We find that 65% of 
subhalos in the range r2oo < r < 2r2oo a-re actually "associated" 
and have thus already been within the host virial radius in the 
past. Roughly one third of subhalos between 2r2oo < r < 3r2oo 
are also physically "associated" with the main halo. The upper and 
lower bounding curves denote the escape velocity for each of the 
five simulated halos. Lower panel: Radial distribution of subhalos. 
Color key is the same as in the upper panel. 



shell gradually settles to equilibrium, approaching a pe- 



4.2. The orbits of associated subhalos 

The discrepancy between the simulation results and 
the naive expectation of the SSIM was pointed out by 
Balogh, Navarro & Morris (2000), and confirmed by sub- 
sequent studies (Mamon et al 2004, Gill et al 2005, Die- 
mand et al 2007) but its physical origin has not yet been 
conclusively pinned down. Associated subhalos found to- 
day beyond their turnaround radius have clearly evolved 
differently from the SSIM prescription, and it is instruc- 
tive to study the way in which the difference arises. 

One possibility is that deviations from spherical sym- 
metry during accretion might be responsible for the out- 
lying associated subhalos. Accretion through the fila- 
mentary structure of the cosmic web surrounding the 
halo, for example, might result in a number of subhalos 
on orbits of large impact parameter that simply "graze" 
the main halo and are therefore not decelerated signifi- 
cantly after their first pericentric approach, as assumed 
in the secondary infall model. These subhalos would lose 
little orbital energy, and should presumably be today 
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Fig. 3. — Orbital trajectories of selected subhalos. Upper panels show the distance to the center of the main progenitor as a function of 
expansion factor. The top-left panel shows the trajectories of 4 subhalos on "extreme" orbits (blue points in Fig.|2j. Note that all of these 
systems gain energy during their first pericentric approach to the main halo. The top-right panel illustrates that interactions occurring 
during the tidal dissociation of bound groups of subhalos are responsible for propelling some satellites onto extreme orbits. At pericentric 
approach, the tidal field of the main halo breaks apart the group, and redistributes each member onto orbits of varying energy. The most 
affected are, on average, the least massive members of the group, some of which are pushed onto orbits with extremely large apocenters. 
The dashed curve shows the growth of the virial radius of the most massive progenitor of the main halo. Bottom panels show the radial 
velocity of the subhalos shown in the upper panels. 



on orbits with apocentric distances of the order of their 
original turnaround radii. According to the analytic cal- 
culation of Mamon et al (2004), systems on such orbits 
may reach distances as far as 2.3r2oo- 

We explore this in Fig. [2 where we show the 
turnaround radius of each associated subhalo versus their 
apocentric distance estimated at z = 0, both normalized 
to the virial radius of the main halo. Subhalos that have 
followed the traditional orbits expected from the SSIM 
should lie to the left of the 1:1 curve in this panel. These, 
indeed, make up the bulk (~ 62%) of the associated pop- 
ulation. 

Note as well that there are a number of subhalos near 
the 1:1 line, whose orbits have not been decelerated since 



accretion into the main halo. These are objects that 
are either on their way to first pericentric passage or, 
as discussed in the above paragraph, that have somehow 
evaded significant braking during accretion. 

More intriguingly. Fig. [2] also shows that there are a 
significant number of subhalos on decidedly unorthodox 
orbits, with apocenters exceeding their SSIM theoretical 
"maximum"; i.e., Tapo > 7'ta- Indeed, 38% of associ- 
ated subhalos are on such orbits, and about ~ 1% are 
on orbits so extreme that their apocenters exceed their 
original turnaround distance by more than a factor of 
~ 2.5 (the latter are highlighted in blue in Fig. [Hif, in 
addition, Tapo > 2r2oo)- The large fraction of systems in 
such peculiar orbits, where the subhalo has gained orbital 
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Fig. 4. — Black dots show the position, at z = 0, of particles 
belonging, at accretion time, to the most massive subhalo of the 
group shown in the right-hand panels of Fig. |3] This subhalo has 
been fully disrupted in the potential of the main halo. Large cir- 
cles show the position of the center of mass of the other (surviving) 
subhalos in the group. Curves show the evolution of each of these 
subgroups since accretion. Note how the surviving subhalos align 
themselves with tidal streams stripped from the main subhalo dur- 
ing the disruption process. The "ejection" of subhalos is thus due 
to the same mechanism that leads to the formation of outgoing 
tidal tails in a merger event and should occur naturally during the 
tidal dissociation of any bound group of subhalos. 



energy since turnaround, indicates that deviations from 
spherical symmetry play a minor role in pushing subha- 
los beyond ^200, and suggests that another mechanism is 
at work. 

4.3. Subhalo mass dependence of unorthodox orbits 

One clue to the mechanism responsible for pushing 
some subhalos onto highly energetic orbits is the depen- 
dence of the effect on the mass of the subhalo. This is 
illustrated in the bottom panel of Fig. [l] which shows 
that low-mass subhalos are the ones being preferentially 
pushed to the outskirts of the halo. 

Further clues result from inspecting individually the 
trajectories of some of the subhalos on extreme orbits. 
This is shown in the top-left panel of Fig. [31 where we 
show the orbits of a few of the associated subhalos with 
''apo > 2.5 rta- Interestingly, all of these subhalos have 
very low mass at accretion (V^^x ~ O.O8V200) and acquire 
their "boost" in orbital energy during their first pericen- 
tric passage. 

The "wiggles" in their trajectories prior to pericenter 
betray the fact that they actually belong to a bound sys- 
tem of multiple subhalos accreted as a single unit (Sales 
et al 2007a, Li & Helmi 2007). This is shown more clearly 
in the top-right panel of Fig. [31 where we show the tra- 
jectories of 5 subhalos belonging to one such group. The 
mass of the group is concentrated in the most massive 
member (see legends in the figure), which is surrounded 



(prior to accretion) by 4 bound satellites. The group 
contributes about 5% of the main halo's mass at accre- 
tion time, Oacc = (1 + ^acc)~^ ~ 0.65. The group as 
a whole turns around at ata ~ 0.35 and accretes on a 
(fpcr:?'apo) = (1:10) orbit that reaches rpor ~ 0.25 r2oo 
at flpcr ~ 0.69. Adding to this evidence, we find that 
~ 95% of subhalos with rapo > 2 r2oo and Tapo > 2.5 rta 
were each, at accretion, members of an FoF group with 
multiple subhalos. 

During pericentric passage, the group is dissociated by 
the tidal field of the main halo, and its 5 members are 
fiung onto orbits of widely different energy. The most 
massive object (single dot-dashed curve in the right pan- 
els of Fig.[3|) follows a "traditional" orbit, rebounding to 
a second apocenter which is only ~ 30% of its turnaround 
distance. The rest evolve differently; the least massive 
subhalos, in particular, tend to gain energy during the 
disruption of the group and recede to a second apocenter 
well beyond the original turnaround. As anticipated by 
the work of Sales et al (2007b) this is clearly the result 
of energy re-distribution during the tidal dissociation of 
the group. 

The bottom panels in Fig.[3l which show the evolution 
of the radial velocity of each subhalo, confirm this sug- 
gestion. The least massive member of the group is, in 
this case, the least bound as well, judging from its ex- 
cursions about the group's center of mass. This subhalo 
(solid black line in the right panels of Fig. [3]) happens to 
be approaching the group's orbital pericenter at about 
the same time as when the group as a whole approaches 
the pericenter of its orbit. This coincidence in orbital 
phase allows the subhalo to draw energy from the in- 
teraction; the subhalo is thus propelled onto an orbit 
that will take it beyond three times its turnaround dis- 
tance, or ^ 6r2oo- Although technically still bound, for 
all practical purposes this subhalo has been physically 
ejected from the system and might be easily confused for 
a system that has evolved in isolation. 

There are similarities between this ejection process and 
the findings of early N-body simulations, which showed 
that a small but sizable fraction of particles are generally 
ejected during the collapse of a dynamically cold N-body 
system (see, e.g., van Albada 1982). The latter occur 
as small inhomogeneities are amplified by the collapse, 
allowing for substantial energy re-distribution between 
particles as the inhomogeneities are erased during the 
virialization of the system. 

In a similar manner, the tidal dissociation of bound 
groups of subhalos leads to the ejection of some of the 
group members. The main difference is that, in this case, 
no major fiuctuations in the gravitational potential of 
the main system occur. Indeed, in the case shown in the 
right-hand panels of Fig.[3]the main halo adds only ^ 5% 
of its current mass and its potential changes little in the 
process. 

A more intuitive illustration is perhaps provided by 
Fig. [H where we show, in the r-VJ-ad plane and at z = 0, 
the location of the same accreted group of subhalos. 
Black dots indicate particles beloging to the main sub- 
halo at the time of accretion. Large circles mark the 
location of the center of mass of each surviving subhalo, 
and the curves delineate their past evolution in the r- 
Vtad plane. The three outermost subhalos track closely a 
stream of particles formerly belonging to the main sub- 
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Fig. 5. — The ratio of apocentric radius (estimated at 2 = 0) 
to turnaround distance as a function of the peak circular veloc- 
ity, Vmax, of a subhalo. Two estimates of Vmax are used for each 
subhalo, one measured at accretion time and another at z = 0. 
Symbols correspond to the median of the distribution, shaded ar- 
eas encompass 25% of the distribution around the median, and the 
extremes of the error bars correspond to the 25th and 75th cen- 
tiles. Note that only fairly massive associated subhalos are today 
on orbits substantially more bound than when they turned around. 
The median apocentric radius of low-mass subhalos is of order of 
the virial radius, indicating that about half of all associated sub- 
halos spend a substantial fraction of their orbital period outside 
''200 • Note that the effect depends only weakly on Vmax below a 
certain threshold; this presumably indicates that, below a certain 
mass limit, subhalos behave like test particles in the potential of 
the main halo. 



halo: the "ejected" subhalos are clearly part of the out- 
going "tidal tail" stripped from the system during first 
approach. The origin of subhalos on extreme orbits is 
thus the same as that of particles at the tip of the outgo- 
ing tidal tails during a merger, and it should therefore be 
a common occurrence during the accretion of any bound 
group of subhalos. 

It is also important to point out that not all low-mass 
subhalos are affected equally. For example, despite being 
of comparable mass to the ejected object, one of the low 
mass members of the group ends up on an orbit almost as 
tightly bound as the main subhalo (red triple dot-dashed 
curve in Fig. [3]). This shows that the orbital fate of a 
subhalo is mainly determined by its orbital phase within 
the accreting group at the time of accretion. Depending 
on this, subhalos may either lose or gain orbital energy 
during the interaction. 

Low mass halos are, however, the ones preferentially 
"ejected" or placed on high-energy orbits through this 
process (see Fig. [1] and Fig. [6]). This is because low- mass 
members of accreting groups will have orbits of larger 
amplitude about its center of mass, enhancing the prob- 
ability of capturing orbital energy when its orbit within 
the group is in phase with the orbit of the group within 
the main halo. In turn, this enhances the survival prob- 



ability of low mass systems by placing them on orbits 
where they spend extended periods in the periphery of 
the main halo, outside the region where tidal fields may 
effectively strip and disrupt them. 

The combination of these two effects (energy gain and 
enhanced survival likelihood) leads to a strong mass de- 
pendence on the orbital properties of associated subhalos 
at z = 0. This is illustrated in Fig.O where we show the 
ratio of apocenter (estimated at z = 0) to turnaround 
distance as a function of subhalo peak circular velocity. 
This figure shows clearly that the most massive subha- 
los are found today in orbits with apocentric distances 
much smaller than their turnaround: halos with V^^x ~ 
0.4 V200 (which corresponds to roughly M^^^ ~ 0.1 M200) 
have median apocenters of order half their turnaround 
distance. On the other hand, the median apocenter of 
associated subhalos with V^^^!^ ^ O.IV200 is of the order 
of the turnaround radius. 

Note as well that the Knax dependence is quite pro- 
nounced at the high-mass end but rather weak for low- 
mass subhalos. This presumably reflects the fact that, 
once a subhalo is small enough, it behaves more or less 
like a test particle in the potential of the main system. 

Finally, note that the mass dependence is less pro- 
nounced when the present-day subhalo VJnax is used. This 
is because tidal stripping has a more pronounced effect 
on systems that orbit nearer the center of the main halo. 
The more massive the subhalo at accretion the closer to 
the center it is drawn and the more substantial its mass 
loss, weakening the mass-dependent bias illustrated in 
Fig. O We shall see below that the mass dependence 
becomes even weaker when expressed in terms of the 
present-day subhalo mass. 

4.4. Subhalo spatial distribution 

The number density profile of all associated subhalos 
is shown by the solid (black) curve in Fig. [6l The pro- 
file may be approximated rather accurately by the same 
empirical formula introduced by Navarro et al (2004) to 
describe the mass profile of CDM halos. This profile is 
characterized by a power-law dependence on radius of 
the logarithmic slope of the density, dlogp/dlogr cx r", 
which implies a density profile of the form, 

ln(n(r)/7i_2) = -(2/a)[(r/r„2)" - 1]- (1) 

This density law was first introduced by Einasto (1965), 
who used it to describe the distribution of old stars 
within the Milky Way. For convenience, we will refer 
to it as the Einasto profile. The scaling parameters n„2 
and r_2 may also be expressed in terms of the central 
value of the density, uq — n(r = 0) = exp(2/Q!) n_2, and 
of the radius containing half of the associated subhalos. 

We list in Table [5] the parameters obtained by fitting 
eq. [l]to the subhalo number density profiles. (Note that 
the units used for uq are arbitrary, but they are consis- 
tent, in a relative sense, for the various subhalo popula- 
tions.) As discussed by Navarro et al (2004), Merritt et 
al (2005, 2006), and more recently by Gao et al (2007), 
ACDM halo density profiles are well described by auM in 
the range ~ 0.15 — 0.3. This is in sharp contrast with the 
much larger values obtained for the subhalo number den- 
sity profile (ofsub ^ 1-0; i.e., the 3D radial distribution 
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Fig. 6. — Number density profile of associated subhalos, after stacking the results of all 5 simulations in our series. The black solid 
symbols show the result for all subhalos, the other symbols correspond to various subsamples obtained after splitting by either V^^'^ (left 
panel) or by subhalo mass at 2 = (right-hand panel). Details on the velocity and mass range for each subsample are given in the legend. 
Solid lines through each curve correspond to the best fits obtained with eq. [T] The parameters of each fit are listed in Table [2] Lines 
without symbols show the dark matter density profile. Note that the spatial distribution of subhalos depends sensitively on subhalo mass 
when measured by V^^xi but that, in agreement with prior work, the mass bias essentially disappears when adopting Alsub to split the 
sample. See text for further discussion. 



of subhalos is approximately "exponential"), and quanti- 
fies the well-established spatial bias between the subhalo 
population and the dark matter mass profile of the main 
halo. The larger values of a characterizing the subhalo 
spatial distribution imply a large nearly constant den- 
sity "core" in their profile, in contrast with the "cuspy" 
density profile of the dark halo, shown as a solid line 
(without symbols) in Fig. [S] 

The left panel in Fig. [H] shows that the subhalo spa- 
tial distribution depends sensitively on subhalo mass, as 
measured by the peak circular velocity at accretion, V;^™ 
(see also, Nagai & Kravtsov (2005), Faltenbacher & r)ie- 
mand (2006), Kuhlen et al (2007)). The various colored 
profiles in this panel correspond to splitting the sample 
of subhalos in four groups, according to the value of Vj^ax 
(normalized to V200 , the virial velocity of the main halo at 
z = 0). The concentration increases systematically with 
Knax! for example, half of the ~ 150 (surviving) subha- 
los with l^ax > 0.17 V200 are contained within ^ 0.7 r2oo 
at 2; = 0. The corresponding radius for subhalos with 
0.04 < V;fiax/^200 < 0.05 is l.lr2oo (see details in 
Table 

Interestingly, the mass dependence of the subhalo 
number density profile essentially disappears when the 
present-day subhalo mass is used to split the sample. 
This is illustrated in the right-hand panels of Fig. [51 
which shows that the shape of the density profile of sub- 
halos differing by up to two decades in mass is basically 
the same. This is in agreement with the earlier results of 



Gao et al (2004) and Diemand et al (2004), but indicates 
that the apparent mass-independence of the subhalo spa- 
tial distribution is not the result of efficient mixing within 
the main halo, but rather a somewhat fortuitous result 
of the cancellation of the prevailing trend by dynamical 
friction and tidal stripping. 

It is conceivable that numerical artifact may also help 
to erase the dependence of nsnhif) on present-day sub- 
halo mass. Indeed, SUBFIND (like every subhalo finder) 
will tend to assign masses to subhalos which depend 
slightly, but systematically, on their location within the 
main halo. The mass of subhalos near the center is more 
likely to be underestimated, and some subhalos may, in- 
deed, even be missed altogether if close enough to the 
central cusp. Splitting the sample by Vt^^x rninimizes 
such effects and allows for the subhalo mass bias to be 
properly established. 

4.5. Velocity anisotropics 

The mass dependence of the subhalo spatial distribu- 
tion discussed in the previous subsection is significant, 
but not very large, and thus is less clearly apparent in 
their kinematics, as shown in Fig. [7l The top panels 
of this figure show the radial velocity dispersion profile, 
Ur = (u^)^/^, computed in spherical shells for the same 
subsamples discussed in Fig. [Sj The bottom panels show 
the anisotropy profile, defined as /3 = 1 — (cr^ -I- cr^)/2cr^. 
The mean values of the velocity dispersion for each com- 
ponent are listed in Tabled 
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Fig. 7. — Radial velocity dispersion and anisotropy profiles for dark matter (thin black lines) and associated subhalos (thick colored 
lines). Symbols are described in the legend and are the same as in Fig. |6] Note that the mass-dependent bias shown in Fig. |6]is also 
reflected in the subhalo kinematics: low mass subhalos tend to have higher velocity dispersions than their high-mass counterparts. This 
bias is clearer when measuring subhalo mass by the peak circular velocity at accretion time, V^'^, rather than by the self-bound mass at 
2 = 0, Mguij. Note as well that subhalos tend to be on orbits less radially biased than the dark matter, especially near the center. This 
is presumably because subhalos on tangentially-biased orbits avoid the innermost regions of the main halo, thus enhancing their survival 
probability. 



The solid lines without symbols in Fig.[7|correspond to 
dark matter particles of the main halo, randomly sam- 
pled in order to match the total number of subhalos. As 
expected, the dark matter velocity distribution is mildly 
anisotropic, with a radial bias that increases outward and 
reaches a maximum near the virial radius. 

The radial velocity dispersion profile of the subhalo 



population follows closely that of the dark matter, al- 
though as a whole, the subhalo population is kinemati- 
cally biased relative to the dark halo. The effect, how- 
ever, is barely detectable; we find a^^^/a^^ ^ 0.98. Our 
results thus confirm the earlier conclusions of Ghigna et 
al (1998), Gao et al (2004), Diemand et al (2004) about 
the presence of a slight kinematic bias between subhalos 
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and dark matter. 

Unlike the conclusions of Diemand et al, however, we 
find a significant discrepancy between the anisotropy pro- 
files of the subhalo population and of the dark halo. As 
shown in the lower panels of Fig. [71 subhalos are on orbits 
less dominated by radial motions than the dark matter 
and, indeed, have a pronounced tangential bias near the 
center (i.e., for r^0.3r2oo). With hindsight, this is not 
entirely unexpected, since subhalos nearer the center are 
more likely to survive if they are on tangentially-biased 
orbits that keep them away from the innermost regions 
of the halo, where tidal effects are strongest. 

4.6. Subhalo mass function 

The large number of associated subhalos on high- 
energy orbits discussed above imply that subhalos within 
the virial radius are just a fraction of all subhalos phys- 
ically influenced by the main halo. This is illustrated 
quantitatively in Fig. [8l where we show the cumulative 
peak circular velocity and mass functions of subhalos 
identified within r2oo- The thin red lines in this figure 
correspond to subhalos identified within r2oo; black to 
the full sample of associated subhalos. Thick lines show 
the average results for the 5 simulated halos considered 
here. The residuals shown in the small panels are com- 
puted relative to the average for the associated subhalo 
population, and show that, on average, the total num- 
ber of associated subhalos exceed those within r2oo by a 
factor of ~ 2. 

Fig. [5] illustrates a number of interesting results. One 
is that, at the low mass end, the shape of the subhalo 
mass and velocity function is insensitive to the radius 
adopted for selection. Indeed, there is no obvious sys- 
tematic trend with Kiax or M^ub for Vmax/^Sax* ^ 0.2. 
Below certain threshold, low mass subhalos behave as 
"test particles" in the potential of the main halo and 
their radial distribution becomes independent of mass. 
This implies that attempts to determine the asymptotic 
slope of the subhalo mass function are unlikely to be 
compromised by selecting for analysis only halos within 
the virial radius, as is traditionally done. 

On the other hand, the subhalo mass function 
shape is substantially affected at the opposite end; 
although about half of all associated subhalos with 
Vmax ^ 0.15 V^max missing from within r2oo, this frac- 
tion decHnes to one quarter for V^ax ~ 0.28 and 
to zero for Vi„ax > 0.31 Vjjjax ■ a result, in that mass 
range, the mass function of subhalos identified within 
^200 is shallower than that of associated systems. This 
should have interesting consequences for semianalytic 
modeling of the luminosity function in galaxy groups 
and clusters, which traditionally assume that all accreted 
subhalos remain within the virial radius of the main sys- 
tem. 

5. SUMMARY AND DISCUSSION 

We have used a suite of cosmological N-body simula- 
tions to study the orbital properties of substructure halos 
(subhalos) in galaxy-sized cold dark matter halos. We ex- 
tend prior work on the subject by considering the whole 
population of associated subhalos, defined as those that 
(i) survive as self-bound entities to z = 0, and (ii) have 
at some time in the past been within the virial radius of 



the halo's main progenitor. Our main findings may be 
summarized as follows. 

• The population of associated subhalos extends well 
beyond three times the virial radius, ^200, and con- 
tains a number of objects on extreme orbits, includ- 
ing a few with velocities approaching the nominal 
escape speed from the system. These are typically 
the low-mass members of accreted groups which 
are propelled onto high energy orbits during the 
tidal dissociation of the group in the potential of 
the main halo. 

• The net result of this effect is to push low-mass 
subhalos to the periphery of the system, creating 
a well-defined mass-dependent bias in the spatial 
distribution of associated subhalos. For example, 
only about ~ 29% of subhalos which, at accretion 
time, had peak circular velocities of order 3% of 
the present-day virial velocity (Vj^^x ~ 0.03V2oo)j 
are found today within r2oo- This fraction climbs 
to ~ 61% and to - 78% for subhalos with Vjf^^^ - 
0.1 V200 and ~ 0.3 V200, respectively. 

• The strength of the bias is much weaker when ex- 
pressed in terms of the subhalo present-day mass, 
due to the increased effect of dynamical friction and 
tidal stripping on the most massive subsystems. 

• The spatial distribution, kinematics, and velocity 
anisotropy of the subhalo population are distinct 
from the properties of the dark matter. Subhalos 
are less centrally concentrated, have a mild velocity 
bias, and are, near the center, on more tangential 
orbits than the dark matter. 

The unorthodox orbits of substructure halos that re- 
sult from the complex history of accretion in hierarchical 
formation scenarios have a number of interesting impli- 
cations for theoretical and observational studies of sub- 
structure and of the general halo population. 

One implication is that subhalos identified within the 
virial radius represent a rather incomplete census of the 
substructures physically related to (and affected by) a 
massive halo. This affects, for example, the interpre- 
tation of galaxy properties in the periphery of galaxy 
clusters, and confirms earlier suggestions that evolution- 
ary effects normally associated with passage through the 
innermost regions of a massive halo, such as tidal trun- 
cation or ram-pressure stripping, should be detectable 
well outside the traditional virial boundaries of a group 
or cluster (Balogh, Navarro & Norris 2000). 

Furthermore, associated subhalos pushed well outside 
the virial radius of their main halo might be erroneously 
identified as separate, isolated structures in studies that 
do not follow in detail the orbital trajectories of each sys- 
tem. This effect would be most prevalent at low masses, 
and it is likely to have a significant effect on the inter- 
nal properties of halos in the vicinity of massive sys- 
tems. We expect, for example, halos in the periphery of 
groups/clusters to show evidence of truncation and strip- 
ping, such as higher concentrations and/or sharp cutoffs 
in their outer mass profiles. 

The same effect may also introduce a substantial envi- 
ronmental dependence in the formation-time dependence 
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Fig. 8. — Associated subhalo mass (Afgub) and peak circular velocity (Vmax) cumulative distributions (both quantities measured at 
2 = 0). Black lines correspond to all associated subhalos; red lines to subhalos identified within r200- Thick lines in each panel denote the 
average of our 5 simulations. Note that the number of associated subhalos exceeds by about a factor of ~ 2 the number of subhalos found 
within r200- The residuals are computed relative to the subhalo population within the virial radius. 



of halo clustering reported in recent studies (Gao et al 
2005; Zhu et al 2006; Jing et al 2007; see also Diemand et 
al 2007b). In particular, at fixed mass, early collapsing 
halos might be more clustered because they are physi- 
cally associated with a more massive system from which 
they were expelled. 

A proposal along these lines has recently been ad- 
vanced by Wang, Mo & Jing (2007) (see also Hahn et al 
2008), who argue that such environmental effects might 
be fully responsible for the age-dependence of halo clus- 
tering. Our physical interpretation, however, differs in 
detail from theirs. Whereas Wang et al argue for the 
suppression of mass accretion onto "old" halos by "heat- 
ing by large-scale tidal fields" as responsible for their 
enhanced clustering, our results suggest that the real cul- 
prit is the orbital energy gain associated with the tidal 
dissociation of bound groups of subhalos, which allows 
"old" low-mass halos to evade merging and to survive in 
the vicinity of massive systems until the present. 

A further implication of our results concern the spa- 
tial bias of the most massive substructures discussed in 
S. 14.41 If, for example, luminous substructures in the Lo- 
cal Group trace the most massive associated subhalos at 
the time of accretion, they may actually be significantly 
more concentrated and kinematically biased relative to 
the dark matter, a result that ought to be taken into ac- 
count when using satellite dynamics to place constraints 
on the mass of the halos of the Milky Way and M31. 

Finally, as already pointed out by Sales et al (2007a, b) , 
gravitational interactions during accretion may also be 



responsible for the presence of dynamical outliers in the 
Local Group, such as Leo I and And XII. Further work is 
needed to assess whether the exceptional orbits of such 
systems could indeed have originated in the tidal disso- 
ciation of groups recently accreted into the Local Group. 
Since the latest proper motion studies of the Magellanic 
Clouds seem to suggest that the Clouds are on their first 
pericentric passage (Kallivayalil et al 2006; Piatek et al 
2007), this is a possibility to consider seriously when try- 
ing to puzzle out the significance of the motion of the 
satellites of the Local Group. 



We thank Simon White and Vincent Eke for useful 
discussions, and an anonymous referee for a constructive 
report. ADL would like to thank Jorge Penarrubia and 
Scott Chapman for many useful discussion which have 
improved this work. The simulations reported here were 
run on the Llaima Cluster at the University of Victoria, 
and on the McKenzie Cluster at the Canadian Institute 
for Theoretical Astrophysics. This work has been sup- 
ported by various grants to JFN and a post-graduate 
scholarship to ADL from Canada's NSERC. AH grate- 
fully acknowledges financial support from NOVA and 
NWO. 



REFERENCES 



Balough M.L., Navarro J.F., Morris S.L., 2000, ApJ, 540, 113 



Bertschinger E., 1985, ApJS, 58, 39 



12 



Ludlow et al. 



TABLE 1 

Properties of simulated halos used in this study. 
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TABLE 2 

Dynamical properties of subhalo population 



Sample 
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-5.1 <Y < -4.6 
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